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Abstract. Assuming that inflation is succeeded by a phase of matter domination, which 
corresponds to a low temperature of reheating T r < 10 9 GeV, we evaluate the spectra of 
gravitational waves induced in the post-inflationary universe. We work with models of hilltop- 
inflation with an enhanced primordial scalar spectrum on small scales, which can potentially 
lead to the formation of primordial black holes. We find that a lower reheat temperature 
leads to the production of gravitational waves with energy densities within the ranges of both 
space and earth based gravitational wave detectors. 
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1 Introduction 

Induced gravitational waves are produced as a result of the interaction between scalar per- 
turbations at second order in the post-inflationary universe. The amplitude of their spectra 
is dependent on the square of the primordial scalar spectrum [1—5], and a relatively large in- 
duced gravitational wave spectrum is expected from the generation of Primordial Black Holes 
(PBHs) [6-9] . In a previous paper Ref. [10] we evaluated the spectra of induced gravitational 
waves generated during a radiation dominated era from the hilltop-type and running mass 
models, which have been shown to be the only models which can lead to Primordial Black 
Holes [11, 12]. We showed that these models lead to an induced gravitational wave signature 
within the sensitivity ranges of planned gravitational wave detectors DECIGO and BBO [13- 
15]. We also found that the running mass model predicted spectra within the sensitivity of 
eLISA [16, 17]under the proviso that inflation is terminated early, with the intriguing factor 
that if we could motivate N <C 40 we would get a detectable signature of PBHs with a mass 
compatible with Dark Matter. 

In this paper, we assume that the universe undergoes a phase of early matter domination 
[18-21] , which lowers the reheat temperature as well as the number of allowed e— folds of 

1 other scenarios are also possible [18, 19, 22, 23] 
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inflation. The source term during the matter phase is constant, and as a result of the 
absence of pressure, the density contrast grows. This may result in perturbations entering 
the non-linear regime and decoupling from the Hubble flow. In our analysis we only assume 
a linear evolution of perturbations, and we account for the non- linear evolution by cutting 
off our analysis at some critical scale. To calculate this critical scale we begin by stating that 
perturbations in the early matter phase evolve as 

^ = -^-V 2 <3? (1.1) 

where p m is the energy density of matter, % is the conformal Hubble parameter and $ is 
the gravitational potential. The evolution of the perturbations is therefore linear until the 
density contrast becomes of order unity which occurs at the scale [24] : 

k NL ~ V^ l, % (1.2) 

where k r is the scale which re-enters the horizon at the time of reheating, is the primordial 
spectrum and k^L is the critical scale at which we terminate our calculation. 

This paper is organised as follows, in section 2 we review the parameters of inflation, in 
section 3 we present the thermal history of the universe, relating the temperature of reheating 
to the relevant scale of reheating, in section 4 we caclulate the bounds on the primordial 
spectrum from PBHs, in section 5 we review the spectrum of induced gravitational waves 
produced during the early matter phase, in section 6 we review the models of inflation that 
can lead to a detectable limit of induced gravitational waves and penultimately in section 7 
we present the results with the final discussion presented in section 8. 

The following conventions are utilised in this paper: r refers to conformal time and 
is related to proper time t as dr = dt/a, a is the scale factor, and the conformal Hubble 
parameter % is related to the Hubble parameter H = a/ a as % = aH. Scales are denoted 
by k, are given in units of inverse megaparsec Mpc -1 and are related to physical frequency 
/ as / = ck/(2air) where c is the speed of light. We assume a radiation dominated universe 
at the time of the formation of the gravitational waves, in which case we have a = ao(T/ro), 
T~L = t , and the scale at re-entry is k = r . 

2 Inflationary Parameters 

Models of inflation can be parametrised by the slow roll parameters [25] : 



s V 

rj = m-pi 



2 



V 

where V is the potential, and derivatives are with respect to the inflaton field (p. These are 
related to the observational parameters, the spectral index n s , the running of the spectral 
index n' s and the scalar spectrum as: 
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n s = 1 + 2r] — 6e 

ri s = 16er] - 24e 2 - 2f 



We use a time re-parametrisation, N = In ( ~J, where the subscripts e and * denote 
the end of inflation and the time of horizon exit respectively. This is related to the potential 
in the slow roll limit as: 

N ~ m pi 2 /""* ^ (2.3) 



and to the scale at horizon exit as [25]: 



N(k ) - N(k) = In (°^\ . (2.4) 

where ko = 0.002Mpc _1 is the pivot scale, and in this paper we effectively take N(ko) = 0. 

We use the latest data release from the WMAP mission [26, 27], for the WMAP data 
combined with BAO and HO data with a null tensor prior. Throughout this paper we take 
n s = 0.96 and n' s < 0.0062. 

3 The temperature of reheating 

The number of e— folds can be related to the temperature of reheating T r as [28, 29] 

where we have taken the energy scale of inflation to be the SUSY GUT scale ~ 10 16 GeV. 
Assuming SUSY means that N max = 56, otherwise we can push this estimate up to N max ~ 
60. In this work we are only interested in reducing the number of e— folds via the inclusion of 
an early matter phase. The thermal history of the universe can support a reheat temperature 
down to lMeV, as this is the temperature below which Neutrinos fail to thermalise and affect 
big-bang nucleosynthesis [30-33] and hence > 37. To relate T r to the scale k r we assume 
the conservation of entropy which gives 



- = ^{g„(n)/g„(Tr)} 1,: ' 



(3.2) 



where g* s is the number of degrees of freedom, Mq is the gravitational scale (Mq = M p / ^/8^T ~ 
2.4 x 10 18 GeV ) and since k = aH/a^ we get the scale which re-enters the horizon at the end 
of reheating k r : 



- 3 - 



4 The spectrum of primordial black holes 



If the primordial spectrum of perturbations towards the end of inflation is large enough, 
i.e if the density contrast exceeds 5 ~ 1/3, perturbations can collapse to form primordial 
black holes. Based on this, constraints can be placed on the spectrum based on astrophysical 
phenomena [34, 35]. In our previous paper, we numerically converted the mass fraction of the 
PBHs into a power spectrum. To perform this analysis we assumed a gaussian distributed 
energy perturbation 2 and a very large reheat temperature T r 3> 10 10 GeV. In this paper, we 
update this calculation for lower reheat temperatures. 

4.1 Mbh(T) and fc(M BH ) 

The comoving wavenumber corresponding to the Hubble radius at temperature T is 

^ -,^16^-1 ( 9*s \V6 ( T 



k = — ~ 1.71 x 10 ib Mpc-M^M 0ntr<l (4.1) 

c y V 106.75/ V lo9GeV AB/ 

On the other hand, the mass of PBHs produced at temperature T is 

M 47r P 3 ^ -1/2^-2 

= 7 y ffl = 4tt 3 / 2 g 3 / 2 1 9 * s T 

« 0.916 x 1O- 2O M ( V ) i^-Y 112 ( q T V 2 , (4.2) 
\(l/V3) 3 ) V106.75/ \W 9 GeVj ' V ; 

where 7 is a numerical factor and Mq = 1.989 x 10 33 g . Eliminating T, we find 

3/4 1/3™ 1/4 
7T ' 5*50 70 5*s 



fc 



45 1 /4G 1 /4 1/3 



2*T y^W7 



1.71 x lO^Mpc" 1 



5*s \-Via / 7 V 7 V M BH " 



(4.3) 



106.75/ V(l/v / 3) 3 y V - 916 x 10-20 M e 



In our numerical calculation we adopt 7 = 1 and g* s = 106.75 for whatever values of 
T , which implies 

M BH « 0.946 x 10 28 g ( 10 ^) 2 , * « 1-67 x 10 9 Mpc" 1 (^g|) " . (4.4) 

These are only precise for Mbh % 10 28 g but the error will not be very large even for larger 
PBHs. 

4.2 T r sets a cut-off 

Our setup is such that the universe after inflation is once dominated by an oscillating scalar 
field and then reheated to the temperature T r . We neglect any PBHs produced before 
reheating (i.e., during the early matter domination) and even those produced after reheating 
if it happens on sub-horizon scales. This gives us conservative PBH constraints. In Ref. [43], 
the authors discussed PBH formation in matter-dominant universe. Assuming a spherical 
collapse of a dense region into a PBH, they obtain the result that more PBHs tend to be 

2 see Ref. [36-42] for alternative scenarios 
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Figure 1. The PBH bound on the primordial spectrum for various cut-off masses, as defined in the 
figure legend. 

produced. This would be reasonable since the pressure P = 0. However, they did not consider 
any non spherical effects, which cause the non-spherical morphologies to evolve during the 
collapsing phase, which would have prevented a further collapse. To get a rigorous bound 
on the PBH spectrum for the early matter phase, we need numerical simulations to obtain 
the correct criterion for matter collapse in matter domination analogous to the radiation 
domination case that 1/3 < 5 < 1 for matter to collapse into a PBH. 

In fact, in this treatment there arises a cut-off PBH mass determined by the reheat 
temperature: 

M BH ,co « 0.946 x 10 32 g (j^—) , (4-5) 
which corresponds to a cut-off wavenumber 



k co ps 1.71 x lf/Mpc- 1 



lGeV 



(4.6) 



The scalar perturbations with k > k co (T r ) may not constrained by PBHs, see Fig. 1. 



5 The spectra of induced gravitational waves 

We follow the analysis of Ref. [4] with much of the notation of Ref. [5]. The spectrum of 
induced gravitational waves generated during an early matter phase was first calculated in 
Ref. [24] and is given as 



V h (k) 



2 /40 N 



\2 poo r\v+l\ v 2 

) dv dy (l-^fv c (ky)V c (kv)I 2 MDS (5.1) 

/ Jo J\v-i\ y 



where x = kr, k is the scale, r is conformal time, v = k/k the ratio of the incoming to 
outgoing scales, y = \fl + v 2 — 2vn, \x is the cosine of the angle between the incoming and 
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outgoing scales, is the spectrum of primordial fluctuations, in this case generated during 
inflation, and Imps is the time integral given as 

3x cos(x) — 3 sin(x) + x 3 
Imds = • (5.2) 

We have taken the lower limit to be x = 0, and terminated the integration at x = kr r , where 
T r is the conformal time at the end of the reheating era. 

The spectrum for induced gravitational waves during radiation domination is given by 

[4,5] 

k 2 poo f \v+l\ 2 

Vh{k) = ^ dv dy^{\ - ^fV c {kv)V Q {ky)hh. (5.3) 

x Jo J\v-i\ y 

where the time integrals, I\ and I2, used in this paper are derived in [10] and given in 
appendix A. 

5.1 Analytical Estimate for a Flat Spectrum 

In this section we assume a flat spectrum and set V^(k) = Ar ~ 10~ 9 . Then the spectrum 
can be written as: 

V h {k) = 2{^\lfdvf dy(l - (l) 2 (5-4) 



consider the term Imps /x 3 ^ 2 : 

Imp _ 

x 3/2 " 



3 

1 H — „ (xcos(x) — sin(x)) 

x 6 



Imps (5.5) 



and it is clear that for x 3> 1 Imps approaches a constant. Substituting this into our equation 
for the spectrum we have 

V h (k) = 2 A% JdvJ dy(l-^f(^j 'l MDS (5.6) 

recall that x = kr, where r is the limit of our time integral which we take to be the end of 
the reheating phase r r = 2/k r . Therefore x r is much greater than 1 for most of the scales we 
consider. Hence we can pull Imps out of the integral 

V h (k) = 2 (f ) 2 A%I MDS JdvJ dy(l - /i 2 ) 2 (j) 2 (5.7) 

and the integrals can be performed analytically. We find that the analytical equation is 
compatible with with numerical calculation for a flat spectrum. 
For scales k r < k < k max the spectrum is then 

Vh ^ /40\ 2 / 16 A; 16k max _ 4k/ 8k r 6 \ 
A R {kf ~ V 3 y V35k max 15 A; 15fc 4 105 A: 6 ; { ' } 

where we have taken the upper limit on v, v max = k max /k and the lower limit to be v m i n = 
kmin/k, and we took k m i n = k r where k r is the scale that re-entered the horizon at the end of 
the matter era, we have also taken Imds ~ 1- By only considering scales which re-enter the 
horizon near k ~ k r , those whose amplitudes have grown the most, Eq. (5.8) can be reduced 
to ~ 356(16A; maz /(15A;)) which for k ma x = k^L ~ 141Av and k ~ k r is 10 5 . Taking instead 
kmax = W 3 k r leads to a spectrum maximum of ~ 10 6 as is confirmed in the full numerical 
calculation shown in Fig. 4. 
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5.2 The evolution of the tensor mode 



Denning = ah k , the equation of motion for the tensor modes is given as 

v k + ( k 2 — — ) = aSk (5. 



and can be solved approximately for the full evolution of the universe. Using step and boxcar 
functions, the source term aSk can be written out as 



aS k oc k 2 r 2 [ 9(r r - r) + ^U TrTeq + (—) 9(r - r eq ) ) (5.10) 



4 / \ 4 

4 li T r r e9 T I 

\ ' V ' eg , f 

where we have taken 5 oc r -3 during radiation domination, 6 is the heaviside step function 
and Ii Tr r eq = Q(j — Tr) — Q(t — T eq ) is the boxcar function. The scale factor can be written 
out in a similar fashion 

a oc r 2 Q(r r - r) + TT r U TrTe + ( — ) t 2 9(t - r eq ) (5.11) 

\ T eqJ 

For sub horizon modes, k 3> 1/r, we obtain the solution plotted black in Fig. 2. In this 
scenario, inflation gives way to an early phase of matter domination which ends when r = r r 
and is followed by a phase of radiation domination that is overtaken by matter at r = r eq . 
The source term is at first constant, then when r r < r < T eq it decays at a rate oc r~ 3 and 
constant again for r > r eq . The amplitude of the sub-horizon tensor modes which re-enter 
the horizon during early matter domination is held at a constant until r r , when begins to 
freely propagate and decay at a rate oc a -1 , until it becomes equal to the source term and 
held at a constant value. The superhorizon modes grow until r = r r , are held at a constant 
between r r < r < r eq and grow again for r > r eq . 

5.2.1 The accuracy of the sudden transition approximation 

Throughout this paper we have utilised the sudden transition approximation between an early 
matter phase and radiation. In this section we investigate the effect a smoother turnover has 
on the tensor modes generated during the early matter phase. For this we make the following 
approximations for the scale factor and source term 



a = 2a r 



\ 2 1 

T \ 1 



T, 



(1 + (r/-r r ) n ) 1 A 



s = ITW (5 ' 12) 

where n is an integer. We plot the results for n = 1, n = 8 and the sudden transition 
approximation in Fig. 3. As is clear from the figure, in all cases the tensor modes approach 
the freely oscillating stage during radiation domination, however, the smooth turnover results 
in a smaller amplitude of an order of magnitude. We also note that n = 8 is very close to the 
sudden transition approximation, and that n = 1 is less than an order of magnitude smaller 
that it. This phenomenon requires further investigation. 
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Figure 2. A schematic of the evolution of tensor modes on different scales with respect to the 
logarithm of the scale factor. The solid red line depicts the evolution of the source term between 
epochs. a r is the scale factor at reheating and a eq is the scale factor at radiation-matter equality. 
Grey lines represent super-horizon modes. The black dash-dot line represents the evolution of a mode 
which enters during EMD. The blue dash-dot line represents the amplitude of the mode which enters 
during RD and the green line, which is barely visible at the right hand side of the plot, represents the 
mode which enters during the current epoch (assuming no acceleration). 




Iog 10 (a) 

Figure 3. The solid lines represent the source term, the red is the sudden transition approximation 
Eq. (5.10), the green is the smooth turnover with n = 1, and the black is n = 4. The dashed lines 
represent the tensor modes. As we can see the tensor modes do approach the freely oscillating limit, 
but there is some loss of amplitude with respect to the sudden transition approximation 
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5.3 Transfer Function 

Detectors of gravitational waves will place a bound on the energy density of gravitational 
waves, denned as [44]: 

Pc d In/ 

where / is the frequency, p c is the critical energy density defining the coasting solution of 
the Friedman equation and pew is the energy density of gravitational waves. This is related 
to the primordial spectrum via a transfer function, Q, = t 2 (k,T)Vh(k). Scales that re-enter 
the horizon during a period of early matter domination experience a constant source term, 
and the amplitude of the tensor mode is kept at it's super horizon value. However once 
the universe enters the radiation epoch the source term decays as does the amplitude of the 
tensor modes, as is depicted in Fig. 2. 

In Ref. [5] they show that for scales smaller than some critical scale, which our scales 
of interest are, the transfer function is 

t(k,r)=^ (5.14) 

where at is the scale factor when the scale k enters the horizon. To be precise, it is the scale 
factor at the time when the source term at that scale begins to decay. That is, for scales that 
enter the horizon during early matter domination, at = a r . Our transfer function is then 

t(fc,r) 



a(r) 
a eq a(r) 

= (5.15) 

where subscript eq is that of radiation-matter equality and the relative energy of scalar- 
induced gravitational waves is 



a eqK e q 

2 

V h (5.16) 



k 



(1 + Z e q) \k r 

where z is the redshift. For scales which re-enter the horizon during radiation domination, 
the relative energy of induced gravitational waves is 

V GW = —*—V h {k) . (5.17) 

J- t %eq 
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5.4 Full numeric results for a flat spectrum 

To get the full spectrum of induced gravitational waves for an early matter phase followed 
by a radiation phase we evaluate 

V h (k) = ( [ T F(v,y,n)dTi + / F(v,y,T 2 )dT2) 

F(v, y, n)F(v, y, T 2 )dT 1 dT 2 + 

f rx T rx rx r fx \ 

/ F(v,y,T\)dT X \ F(v,y,T 2 )dr 2 + / F(v,y,T 2 )dr 2 / F(v, y, r^dn J 

\J J x r JO J x r J 

/■x 

+ J F(v,y,T 2 )F(v,y,n)dTidT 2 

.matter + c h .cross + v h .radiation 
« v h .matter + v h .radiation 

(5.18) 

where F(v,y,r) is the integral over v,y in Eq. (5.3), and we drop the cross terms arising 
from < $ matter^ radiation > hi the last line. This approximation is reasonable since the cross 
terms are only of significance at k ~ k r . 

Figure 4 is a depiction of the spectrum of induced gravitational waves arising from a 
flat primordial spectrum of density perturbations; n s = 1. This figure was generated mainly 
for illustrative purposes, and as such we have chosen kj^L to be 10 3 times as large as the 
value calculated using Eq. (1.2). We assume that modes with k > k^L do not experience 
the constant source term. Scales which are still super-horizon at the end of the early-matter 
phase have a spectrum Vh oc A; 3 and therefore become rapidly smaller than the spectrum 
generated by the pure radiation source term. That means that for modes which enter the 
horizon soon after r r we need only consider the convolution of modes with those that enter 
during the radiation era. 

We present the results for a flat spectrum at various reheat temperatures in Fig. 5, 
where we have taken the limits on v to be v m i n = k r /k and v max = k^L/k with the latter 
upper bound accounting for the non-linear cutoff. We could have modified the calculation 
and checked for each v and y that k < k^L and |k — k| < &tvl, however simply modifying 
the limits of v has the same effect. 
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Figure 4. We plot the spectrum of induced gravitational waves for a flat primordial spectrum with 
T r = 10 9 GeV. The spectrum for scales that re-enter the horizon deep in the radiation era have a flat 
spectrum, due to the fact that the source term is a decaying function and thus the modes oscillate 
freely, and is represented by the solid blue line. The solid black line represents the modes that re- 
enter the horizon during the early matter phase (k > fc#). The The red dashed line is the complete 
spectrum, assuming an early matter phase followed by a phase of radiation domination. The spectrum 
for k > fcjvL behaves as Vh oc 1/fc 4 , as Vh oc l/k for k r < k < Hnl, as Vh oc fc 3 for k < k r and 
as Vh ~ constant for k <C k r . One can think of this as follows: modes that re-enter the horizon 
during the radiation phase but with k ~ k r , i.e. near the EMD phase, will interact with modes 
that re-entered during EMD and hence their behaviour/characteristics are modified from the instant 
reheating scenario. We have utilised a simplified analysis, in that Vh(k) = 'Ph matter (k) + r Ph rad (k). 
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Figure 5. We plot the spectrum of induced gravitational waves for a fiat primordial spectrum 
with various reheat temperatures and the cutoff scale k max — fcjVL- The black dashed lines cor- 
respond to the PBH bound assuming (from left pseudo-vertical line to right pseudo-vertical line) 
T r = 1 GeV, fOGeV, 10 2 GeV, f0 3 GeV, f0 4 GeV, 10 5 GeV, f0 6 GeV, f0 7 GeV, f0 8 GeV and 10 9 GeV. The 
cluster of solid lines in the top right corner are the sensitivity ranges of ground based detectors, L1G0 
S5 and S6 [45], and KAGRA [46], while the thick horizontal salmon pink line is the forecast sensi- 
tivity of Advanced LIGO [47, 48] . Also shown is the sensitivity limit of the Square Kilometre Array 
(SKA) [49-52]. The green, red, blue and purple solid lines correspond to taking T r — lMeV, lGeV, 
f 4 GeV, and 10 8 GeV. The blue dashed line is the spectrum for T r = 10 4 GeV without terminating at 
fcjvL • It is interesting to note that even a flat primordial spectrum can lead to a spectrum of induced 
gravitational waves detectable by cross-correlated DECIGO 
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6 The Models of Inflation 



The spectrum of induced gravitational waves is directly proportional to the square of the 
primordial spectrum and is therefore clear that for an enhanced spectrum of induced grav- 
itational waves one needs an enhanced primordial spectrum. Since at the pivot scale the 
spectrum is tightly constrained by CMB data we need to go beyond this and consider models 
which enhance the spectrum on small scales. Phenomenologically, the two models of infla- 
tion that exhibit this property [11, 12] are the running mass model and the hilltop model, 
depicted in Fig. 6. 

6.1 The Hilltop type model 

Identified in Ref. [11] as the phenomenological form necessary for PBH formation this model 
has the potential [53, 54]: 

V = Vo(l + ripifP — rj q ip q ) (6.1) 

where the coupling terms r/ p and r\ q are less than 1 and p < q to get the shape in Fig. 6. 
Certain realisations in super gravity can be found; see for example Refs. [11, 55-60]. 

This model is very compatible with WMAP data, and as such many of the model's 
terms are not ruled out, but in this analysis we add the extra requirement that the model 
is maximised at small scales and yet still remain within the PBH bound. We also take the 
basic number of e— folds to be N = 56. Parameter selection criteria is explained in more 
detail in Ref. [10]. 

6.2 The Running Mass Model 

This model is the basic 4> 2 model with a varying mass term that arises in taking renormalised 
group equations, and is given as [8, 61-68] 

V B 2 , Aip 2 

Vb =1 "T^ + 2(l + ahnW (6 - 2) 

In this case we select parameters which satisfy n s = 0.96, n' s = 0.0039 and n' s = 0.0043, 
which for T r > 10 9 GeV are terminated at N = and N = respectively. 



Figure 6. An illustration of the Hilltop- type and running models. In our scenario scales of cosmo- 
logical interest during the hilltop regime (indicated with red), and the end of inflation occurs once 
the inflaton has reached a flatter region of the potential (indicated with blue) . 
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7 Results 



We plot results for the hilltop model for p = 2 and q = 2.3,3,4, and for the running mass 
models which satisfy n s = 0.96 and n' s < 0.0062 with N ~ 57. These are plotted in Fig- 
ures. 7,8, 9, 10, 11 and 12 for a range of reheat temperatures lGeV < T r < 10 9 GeV 3 . Each 
reheat temperature modifies the maximum allowed number of e— folds N max and therefore 
we have integrated only up to k max = k p i vo te Nmax , except for in the case when k max > k^L 
where we only integrate up to k^L- In the final figure, Fig. 12 we have plotted the results 
of the hilltop and running mass models for a reheat temperature of T r = 10 6 GeV. 

In our previous paper Ref. [10], we calculated the spectrum of Induced Gravitational 
Waves for the running mass models with large running 0.0067 < n' s < 0.012 which is no longer 
supported by the latest WMAP release [26, 27]. On a related note, to motivate N < 37 by 
modifying the reheat temperature would require T r < lMeV which is unsupported by theory. 

We also find that for the running mass model satisfying N = 38 e— folds, the corre- 
sponding induced gravitational waves spectra were not within the sensitivity ranges of any 
of the experiments. However, if the k^L cutoff can be relaxed, the spectra may very well be 
within the ranges of SKA and PULSAR. 



3 There is an error in our previous paper, a missing factor of 4 in the / function which means that our 
previous results are 4 2 times smaller than they should be, this is corrected for here. 




Figure 7. The spectra of induced gravitational waves for the Hilltop model with p — 2 and q = 3, 
with k max — kjsiL- The black dashed lines and the cluster of solid lines in the right hand corner are 
defined in Fig. 5. The red line represents the induced gravitational wave spectrum for T r ^> 10 10 GeV 
and N = 55 e— folds. From right to left the black solid lines correspond to the induced gravitational 
wave spectrum for reheat temperatures of T r = 10 8 GeV, T r — 10 5 GeV, and T r = lGeV. 
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Figure 8. The spectra of induced gravitational waves for the Hilltop model p = 2 and q = 4, with 
k m ax = kjsrL- The black dashed lines and the cluster of solid lines in the right hand corner are defined 
in Fig. 5. The red line represents the induced gravitational wave spectrum for T r 3> 10 10 GeV and 
N = 55 e— folds. From right to left, the black solid lines correspond to the induced gravitational wave 
spectrum for T r = 10 8 GeV, T r = 10 5 GeV, and T r = lGeV. For illustration, we include the spectrum 
for a reheat temperature of T r = 10 J GeV and integrated up to the maximum scale instead of the 
non- linear cutoff (dashed blue line). 
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Figure 9. The spectra of induced gravitational waves for the Hilltop model p = 2 and q — 2.3, with 
kmax — fcjVL- The red line represents the induced gravitational wave spectrum T r ^> 10 10 GeV and 
N = 55 e— folds. From right to left, the black solid lines correspond to the induced gravitational wave 
spectrum for T r = 10 8 GeV, T r = 10 5 GeV, and T r = lGeV. The black dashed lines and the cluster of 
solid lines in the right hand corner are defined in Fig. 5. 
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Figure 10. The spectra of induced gravitational waves for the running mass model satisfying n' s — 
0.0039 and N = 57 , with k max = Rnl- From right to left, the black solid lines correspond to the 
induced gravitational wave spectrum T r — lGeV,T r = 1 x 10 5 GeV, and T r — 10 9 GeV. The black 
dashed lines and the cluster of solid lines in the right hand corner are defined in Fig. 5. 
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Figure 11. The spectra of induced gravitational waves for the running mass model satisfying n' s = 
0.0043 and N — 57, with k max = fcjv2> From right to left, the black solid lines correspond to the 
induced gravitational wave spectrum T r — lGeV,T r = 1 x 10 5 GeV, and T r — 10 8 GeV.T The black 
dashed lines and the cluster of solid lines in the right hand corner are defined in Fig. 5. 



- 18 - 




Figure 12. The spectra of induced gravitational waves for all the models plotted in the previous 
figures assuming a reheat temperature of T r = 10 5 GeV, with k max — fcjvx- The solid lines are hilltop 
with p = 2 q = 2.3 (red), p — 2 q — 3 (black) and p = 2 q = 4 (blue); the dashed lines are the running 
mass models with n' s = 0.0043 (blue) and n' s — 0.0039 (red). The black dashed lines and the cluster 
of solid lines in the right hand corner are defined in Fig. 5. 
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8 Discussion 



In this work we assumed a sudden transition between early matter domination and radiation 
domination. In the early universe however, the Hubble time can be of the order of the decay 
rate of matter to radiation. In this case modes re-entering the horizon towards the end of 
matter domination will also experience a decaying source term and the effect shown here 
may be reduced, depending on how long the transition phase lasts, as is shown in Fig. 3. 
Therefore the results presented here are upper bounds on what the induced gravitational wave 
spectrum could be, with the actual spectrum possibly being only an order of magnitude or 
two smaller than what we calculated. 

Therefore our conclusions depend on the condition that the tensor modes generated 
during an early matter phase survive the transition into radiation. Under this proviso, we 
have shown that assuming an early matter phase with T r = 10 5 GeV results in a spectrum 
of induced gravitational waves with energy densities within the range of BBO/DECIGO and 
cross-correlated DECIGO [69] . Since the assumption of an early matter phase truncates the 
PBH bound at a smaller value of A:, we have shown that the running mass model generates 
induced gravitational waves detectable by ground based gravitational wave detectors such as 
LIGO, and KAGRA. This means that the RMM model with a running of n' s = 0.0043 and a 
reheat temperature of 10 8 GeV as well as the model with n' s = 0.0039 and T r = 10 9 GeV can 
be ruled out since LIGO has failed to detect a gravitational signature. 
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A The time integral during the radiation era 

if 4 4 } 

h = ^ ky 3 v 3 | ~ cos(x) «nSi(/? n x) + sin(x) y^(-l) n+1 q n ci(/3 n x) j 

+71 sin(x) + 72 sin(vx) sin(yx) + 73 sin(vx) cos(yx) + 74 cos(vx) sin(yx) + 75 cos(vx) cos(yx) 

(A.l) 

and the second integral is given by: 

h = — 7T7~~3~~3 {cos(x) [-Si^x) + Si(/3 2 x) + Si(/3 3 x) - Si(/3 4 x)] 
+ sin(x) [ci(/3ix) + ci(/?2x) - ci(/3 3 x) - cif^x)]} 

+721 sin(x) + 722 sin(x(v + y)) + 723 sin(x(u - y)) + 724 cos(x(v - y)) + 725 cos(x(v + y)) 

(A.2) 

the coefficients in these two integrals are given in the tables below. The a n coefficients in the 
first integral satisfy Xm=i( — l) nQ; n = an d the a coefficient in the second integral is given 
by a = (v 2 — 1 + y 2 ) 2 . 
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The Si and ci terms are the sine and cosine integrals respectively [70], defined as: 



Sifx) 



CI X 



sin(i) 
x cos(t) 



dt 



dt . 



(A.3) 
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-1 + v + y 




1 + v — y 


/3 4 


— 1 + v — y 



Coefficient 
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sin(wx) sin(yx) 
s'm(vx) cos(yx) 
cos(vx) sin(yx) 
cos(vx) cos(yx) 



sin(x) 
sin(x(v — y)) 
sin(x(v + y)) 
cos(x(v — y)) 
cos{x{v + y)) 
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Expression 
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kx 3 l y 3 v 3 (-2 + x 2 - x 2 y 2 - xV - 2x 2 t>y) 
- jggg (2 - x 2 + s 2 y 2 + xV - 2x 2 i»y) 



Table 1. Table on the left gives the expressions for the coefficients of the arguments of the Cosine 
and Sine integrals in Eqs. (A.l) and (A. 2). The right table gives the expressions for the t\ integral 
Eq. (A.l) (tob block) and for the Ti integral Eq. (A. 2) (bottom block). 
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Table 2. This table gives the expressions of the coefficients of the sine and cosine integrals in 
Eq. (A.l). Each a n coefficient has the same parameters as the others, but the parameters differ in 
their respective signs. The columns to the right of the as give the sign of the parameter denned in 
the column header. For example then we can read off o,\ as — 1 + v 4 + 4i> 3 + 4t> 2 — 3y 4 — 4y 3 — 2y 2 v 2 . 
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